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Multi-model assessment of the factors driving 
stratospheric ozone evolution over the 21 st century 
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45 Abstract 

46 The evolution of stratospheric ozone from 1960 to 2100 is examined in simulations from 

47 fourteen chemistry-climate models. There is general agreement among the models at the 

48 broadest levels, showing column ozone decreasing at all latitudes from 1960 to around 2000, 

49 then increasing at all latitudes over the first half of the 21 st century, and latitudinal variations in 

50 the rate of increase and date of return to historical values. In the second half of the century, 

51 ozone is projected to continue increasing, level off or even decrease depending on the latitude, 

52 resulting in variable dates of return to historical values at latitudes where column ozone has 

53 declined below those levels. Separation into partial column above and below 20 hPa reveals that 

54 these latitudinal differences are almost completely due to differences in the lower stratosphere. 

55 At all latitudes, upper stratospheric ozone increases throughout the 21 st century and returns to 

56 1960 levels before the end of the century, although there is a spread among the models in dates 

57 that ozone returns to historical values. Using multiple linear regression, we find decreasing 

58 halogens and increasing greenhouse gases contribute almost equally to increases in the upper 

59 stratospheric ozone. In the tropical lower stratosphere an increase in tropical upwelling causes a 

60 steady decrease in ozone through the 21 st century, and total column ozone does not return to 

61 1960 levels in all models. In contrast, lower stratospheric and total column ozone in middle and 

62 high latitudes increases during the 21 st century and returns to 1960 levels. 

63 

64 
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64 1. Introduction 

65 

66 The evolution of ozone in the 21 st century is a critical science issue. While changes in 

67 ozone are presently controlled primarily by declines in halogen concentrations, variations in 

68 temperature, circulation, and oxides of nitrogen and hydrogen also affect ozone [WMO 2003; 

69 WMO 2007]. Throughout the stratosphere there are long-term changes in various processes and 

70 the balances among them. Therefore it is difficult to find a single approach to identify the 

71 contributions of different mechanisms affecting ozone levels throughout the stratosphere. 

72 Ozone loss throughout much of the stratosphere was dominated by rising halogen 

73 concentrations until about 1997, the time of peak halogen loading [e.g., Yang et al., 2006; 

74 Shepherd and Jonsson, 2008; Yang et al., 2008]. In the 21 st century, as halogen concentrations 

75 are expected to decrease at about 1/3 of the rate of their increase in the late 20 century, other 

76 factors will likely play a more significant role in the evolution of ozone. Stratospheric cooling 

77 will decrease the rate of gas-phase reactions that destroy ozone, and thereby increase 

78 concentrations of ozone [e.g., Haigh and Pyle , 1979; Brasseur and Hitchman, 1988; Shindell et 

79 al., 1998; Rosenfield et al., 2002], Future increases in N 2 0 and CH 4 increase loss of ozone by 

80 nitrogen and hydrogen catalytic cycles, [e.g., Randeniya et al, 2002; Rosenfield et al., 2002; 

81 Chipperfield and Feng, 2003; Portmann and Solomon, 2007; Ravishankara et al., 2009] but in a 

82 changing climate, increases in N 2 0 do not necessarily cause NO y (or NO x ) to increase 

83 proportionately throughout the stratosphere. Rosenfield and Douglass (1998) showed that NO y 

84 loss rates increase as the stratosphere cools, which will impact future levels of ozone. Several 

85 studies have linked increases in greenhouse gases (GHGs) to changes in stratospheric transport. 
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86 further impacting ozone recovery [Waugh et ai, 2009; Li et ai, 2009; Hegglin and Shepherd, 

87 2009], 

88 Eyring et al. [2005] laid the groundwork for the Chemistry-Climate Model Validation 

89 (CCMVal) activity in which numerous chemistry-climate models (CCMs) were evaluated to 

90 increase our confidence in predicting future stratospheric ozone change. The results of the 

91 CCMVal-1 simulations of past ozone changes conducted for this activity were presented in 

92 Eyring et al. [2006], which evaluated processes important in determining the distribution of 

93 ozone. In a follow on study Eyring et al. [2007] presented the projections of stratospheric ozone 

94 change over the 21st century simulated by these CCMs, and discussed how quantities that can 

95 impact ozone are projected to change. 

96 CCMVal- 1 was the first major coordinated activity where an ensemble of CCMs 

97 performed simulations with similar external forcings to assess ozone evolution [ Eyring et al., 

98 2007], The CCMVal-2 activity [. Eyring et al., 2008] included more models (14 are used here) 

99 and, in addition, more simulations cover the entire period of interest (1960-2100) than in 

100 CCMVal- 1 (3 models). Here, we make use of CCMs that participated in the CCMVal-2 activity 

101 and that contributed projections of stratospheric ozone evolution until the end of the 21 st century. 

102 These simulations are based on observed GHG and halogen concentrations in the past and on one 

1 03 predicted scenario for the future. 

104 Austin et al. [2010a] examined the evolution of total column ozone and compared ozone 

1 05 and Cl y recovery dates in the CCMVal-2 models. Here we (a) contrast the evolution in the upper 

106 and lower stratosphere, and (b) examine the causes of the ozone changes (and differences among 

107 the models). To accomplish this, we apply a variety of methods depending on the physical 

108 and/or chemical processes that dominate ozone evolution in various regions. Section 2 gives an 
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109 overview of the models and model simulations. The results are discussed for extra-polar ozone 

1 10 in Section 3 and polar ozone in Section 4. Variations in ozone recovery by region are discussed 

111 in Section 5 and a summary of the conclusions is given in Section 6. 

112 

113 

114 2.1 Models and Model Simulations 

115 

116 CCMVal-2 model simulations were conducted to improve the understanding of models 

117 through process-oriented evaluation along with discussion and coordinated analysis [SPARC, 

118 2010]. Information about individual model simulations and references for each are summarized 

119 in Table 1. In addition, Morgenstern et al. [2010] present a much more detailed overview of the 

120 models that participated in this activity. Here we consider 1960 to 2100 simulations from the 

121 CCMs whose domains include at least the upper stratosphere, and can then be used to contrast 

122 the ozone evolution in the lower and upper stratosphere. All models in this study have output for 

123 1960-2099, except UMUKC A-MET O which ends in 2083. 

124 REF-B2 uses the Alb GHG scenario from IPCC [2000] and the adjusted Al halogen 

125 scenario from WMO [2007]. The adjusted Al halogen scenario includes the earlier phase out of 

126 hydrochlorofluorocarbons mandated by recent revisions to the Montreal Protocol. For 

127 GEOSCCM a combination of REF-B1 (1960-2000) and REF-B2 (2001-2099) was used. REF- 

128 B1 differs only in that observed Hadley sea surface temperatures (SST) and sea ice data [Rayner 

129 et al., 2003] were used instead of modeled SSTs. All model runs with the exception of CMAM 

130 use prescribed sea surface temperatures and sea ice extent based on fully coupled atmosphere- 

131 ocean runs for the Alb scenario. CMAM includes a fully coupled ocean model in its 
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132 simulations. All models use sophisticated stratospheric chemistry schemes except AMTRAC3, 

133 in which halogens are not explicitly modeled but rather Cl y and Br y are parameterized [Austin 

134 and Wilson, 2010; Morgenstern et al., 2010]. REF-B2 simulations use background non-volcanic 

135 aerosol loading and there are no imposed solar cycle variations. Although several modeling 

136 groups submitted multiple ensemble simulations, only one member from each group was used in 

137 this analysis because the intermodel differences are generally much larger than the ensemble 

138 spreads [Eyring et al., 2007]. The model output used in this study is a zonal mean monthly mean 

139 average. The partial and total column ozone amounts are integrated over the standard CCMVal- 

140 2 pressure levels. 

141 In order to view all 14 models on a single plot, filtering is necessary to remove short-term 

142 (e.g. high frequency) variations. In all cases (except where noted) a 1:2:1 filter is used iteratively 

143 30 times as described in Eyring et al. [2007] to smooth the model output displayed in the figures. 

144 While some differences can be seen in individual models between this filter and the time series 

145 additive-model (TSAM) analysis, they are typically very small and do not impact our 

146 conclusions [Scinocca et al. 2010, in preparation], 

147 

148 2.2 MLR Analysis 

149 

150 One of the primary methods we use to estimate the contribution of different mechanisms 

151 to the simulated changes in ozone is multiple linear regression (MLR). The method is explained 

152 in detail in Oman et al. [2010], but we repeat some of the basic details here. For a given location 

1 53 and time, MLR is applied to determine the coefficients m x such that 

AO 3 (0 = A x j(t) + 


154 


( 1 ) 
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155 where the X f are the different quantities that could influence ozone, the coefficients m x are the 

156 sensitivity of ozone to the quantities X , i.e., m x = d0 3 /dX , and s is the error in the fit. To do 

1 57 this four explanatory variables ( Xj ) are used in ( 1): (i) Cl y + a Br y , (ii) reactive nitrogen (NO y 

158 = N O+N O2+NO3+2 * (N 2 0 5 )+HN0 3 +H02N O2+CION 0 2 +Br0N 0 2 ), (iii) reactive hydrogen (HO x 

1 59 = OH+HO2), and (iv) temperature (T). Each term on the right hand side of equation (1) gives the 

160 “contribution” of the response in ozone due to a change in X . We use a = 5 in the definition as 

161 Daniel et al. [1999] show this is an appropriate value for the upper stratosphere, which is the 

162 chemically driven region for which the MLR analysis is most appropriate. A test of this method 

163 using a = 60 did not noticeably impact any of the individual contributions calculated in this 

164 study. We do use a = 60 in section 4 when examining polar lower stratospheric ozone changes. 

165 There are several limitations with the above linear regression approach as discussed in 

166 Oman et al. [2010] which are summarized below. First, other mechanisms that are not considered 

167 in the regression (e.g., transport) could play a role. Second, significant correlations can exist 

168 between the temporal variations of the quantities, i.e., the quantities are not necessarily 

169 independent. Third, a high correlation between ozone and a quantity does not show causality, as 

170 ozone could be causing the quantity to change, or changes in another quantity could be causing 

171 both ozone and the quantity of interest to change in a correlated way. Temperature and ozone in 

172 the upper stratosphere are an example of this third complication: changes in ozone cause changes 

173 in temperature through changes in short-wave heating [Shepherd and Jonsson, 2008], At the 

174 same time, the local ozone concentration responds to changes in temperature by changing 

175 chemical reaction rates. Also, the relationship between the regression variables and ozone may 

176 not be linear. For example, the C10/Cl y ratio varies as a function of temperature in a highly non- 

177 linear manner. One could use a regression based on the rates of odd-oxygen loss by nitrogen, 
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178 chlorine, and bromine radical abundances as regressor variables, as was done by Yang et ah, 

179 [2006], However, this treatment is beyond the scope of the present analysis, as the information 

180 archived by each model makes use of radical abundances challenging. Because of the above 

1 8 1 limitations, caution must be applied when interpreting the MLR results presented below. 

1 82 UMUKCA-UCAM and AMTRAC3 did not provide HO x so for these models the MLR 

183 consists of three explanatory variables. This analysis has been repeated without HO x for all 

184 models and while there is some decrease in explained variance from the MLR the main findings 

185 and conclusions are not impacted. 

186 

1 87 3. Extra-Polar Ozone 

188 

189 3.1 Total column ozone 

190 

191 Austin et al. [2010a] examined the changes in total column ozone in the CCMVal-2 

192 simulations. Here, we start with total column ozone and then break down the ozone amounts 

193 into two partial columns for the extra-polar region. Figure la-c shows the evolution of total 

194 column ozone amounts over the tropics (25°S-25°N) and midlatitudes of each hemisphere (35- 

195 60°S and 35-60°N), for each of the CCMVal-2 models. The evolution is shown with respect to 

196 1960 levels, and has been smoothed as described above. Also shown are ground-based total 

197 column ozone measurements (solid black curves, updated from Fioletov et ah, [2002]) with 

1 98 respect to the smoothed 1 964 value. 

199 There is qualitative agreement in the evolution of ozone between regions and models at 

200 the broadest level, with ozone decreasing from 1 960 to around 2000, and then increasing in the 
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201 first half of the 21 st century [Austin et al., 2010a]. All models consistently indicate tropical total 

202 column ozone at the end of the 21 st century is less than it was in 1960 (Figure lb). In contrast, 

203 total column ozone over the midlatitudes (Figures la,c) in 2100 exceeds 1960 values in nearly all 

204 the models, with a generally larger increase in Northern Hemisphere (NH) mid-latitudes. 

205 There is a large spread in the magnitude of the peak ozone loss (around 2000) and the 

206 date that column ozone is predicted to return to 1960 levels. For example, in SH mid-latitudes 

207 the simulated decrease in ozone from 1960 to 2000 varies from around 10 DU to over 50 DU, 

208 while the date of return to 1960 levels for the same region varies from around 2030 to after 2100. 

209 In all three latitude bands CNRM-ACM and MRI show the largest loss in total column ozone by 

210 the year 2000, and also tend to have later return dates. UMUKCA-METO has the smallest ozone 

211 loss in the tropics and Southern Hemisphere (SH), and one of the smallest decreases in NH 

212 midlatitudes. The reasons for these model differences are examined in Section 3.2. Also, there 

213 is a large difference in the ozone decline in the SH midlatitudes prior to 1975 found by the 

214 various models (declines ranging from 0 to 20 DU), with no change over this time period 

215 apparent in the ground-based ozone observations. Ground-based total column ozone observations 

216 indicate about a 20 DU decrease by 2000 in the SH midlatitudes and about a 15 DU decrease in 

217 the NH midlatitudes. In the tropics only small total column ozone changes have been observed 

218 and are less than the model-simulated changes with the exception of UMUKCA-METO, which 

219 are similar. 

220 

221 3.2 Partial column ozone 


222 
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223 While the evolution of total column ozone and its impact on surface UV radiation 

224 changes are of great importance, it is useful to look at changes in the partial columns [e.g., Yang 

225 et al., 2006; Li et al., 2009J. This shows the relative role of changes in the upper and lower 

226 portions of the stratosphere in the differences described above. As in Oman et al. [2010] we split 

227 the stratosphere at 20 hPa because the ozone changes at 20 hPa are generally small, and 20 hPa 

228 separates the photochemically controlled ozone region above from the transport and chemically 

229 driven region below. Other studies have used 15 hPa [e.g. Li et al., 2009] or 25 km altitude (~25 

230 hPa) [e.g., Yang et al., 2006] but we prefer 20 hPa, as used in Oman et al., [2010]. 

231 Figures ld,e,f and lg,h,i show the evolution of the partial column ozone amounts for the 

232 upper (20 to 0.1 hPa) and lower (500 to 20 hPa) portions of the column. The evolution in the 

233 upper partial column ozone is similar among the 3 regions, but there are significant differences in 

234 the lower partial column. It is the differences in the lower stratosphere that cause differences in 

235 the evolution of column ozone between regions. In particular, the decrease in tropical column 

236 ozone over the latter half of the 21 st century is due to decreases in the lower stratosphere (Figure 

237 lh), which are larger than increases in upper stratospheric ozone (Figure le). 

238 There are quantitative differences in the upper columns but these are generally smaller 

239 than the differences in the lower columns, especially in midlatitudes (note the different scales for 

240 upper and lower column plots in Figure 1). Peak losses occur around the year 2000, averaging 3 

241 to 7 DU over tropical and midlatitudes for the upper partial column of most models, with MRI 

242 and CNRM-ACM having losses around 10 DU. In the lower partial column there are much 

243 greater differences in the ozone change over this range of latitude bands. Midlatitude ozone 

244 losses around 2000 are largest in the SH, as expected due to the effects of the Antarctic ozone 

245 hole [Atkinson et al., 1989], Ozone losses range from 10-25 DU with less ozone loss in 
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246 UMUKCA-METO and significantly more in MRI and CNRM-ACM. Tropical latitudes show 

247 largest loss of ozone by the end of the 21 st century, with a range of 7 to 18 DU. The tropical 

248 lower partial column is the only region where ozone is consistently forecast to continue to 

249 decline until the end of the 21 st century, due to increased strength of the Brewer-Dobson 

250 circulation in the models that is driven by rising GHG concentrations [Rind et al., 1998; Buchart 

251 and Scaife, 2001; Butchart et al., 2006; Shepherd, 2008; Li et al.. 2009], It should be noted that 

252 models with larger (smaller) total column ozone changes generally have larger (smaller) changes 

253 in both lower and upper partial columns. 

254 

255 3.3 Upper Stratosphere 

256 

257 As discussed above, and shown in Figure ld-f, the evolution of ozone in the upper portion of the 

258 stratosphere is qualitatively similar among most models. Also, the overall behavior of ozone in 

259 the tropics and midlatitudes (both hemisphere) is similar. Nonetheless, there are quantitative 

260 differences among the models. We use the MLR analysis described in Section 2 to examine the 

261 mechanisms causing the ozone changes, and the differences among the models. 

262 We consider first the evolution of ozone at 5 hPa in the tropics (25°S-25°N). Figure 2a 

263 shows the evolution of the ozone mixing ratio while Figure 2b shows the change in ozone with 

264 respect to 1960 levels. There is a large spread in the time-mean values of ozone among the 

265 models (from around 8 ppmv to around 1 1 ppmv), and also some differences in changes relative 

266 to 1960 levels. All models show a peak ozone loss around 2000, with magnitudes of 0.5 and 0.7 

267 ppm for most models. The MRI and CNRM-ACM models show the largest peak ozone loss 

268 (greater than 1.0 ppm), whereas AMTRAC3 and LMDZrepro have the smallest loss (less than 
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269 0.4 ppm). All models also show 5 hPa ozone returning to 1960 levels before the end of the 21 st 

270 century, but there is a large spread in dates when this occurs (from early 2020s to 2060s). 

271 As presented in the Introduction, several mechanisms can cause changes in ozone, 

272 including changes in halogens, temperature, reactive nitrogen and hydrogen, and transport. 

273 Figures 2c,e,g,i show the evolution of Cl y +aBr y , temperature, NO y , and HO x , respectively. The 

274 evolution of Cl y +aBr y and variation among models (Figure 2c) is somewhat similar to that of 

275 ozone (Figure 2a), suggesting that differences in Cl y +aBr y can explain much of the ozone 

276 evolution. However, there is not a simple one-to-one relationship between changes in ozone 

277 and those in Cl y +aBr y . For instance, ozone exceeds 1960s values in the later part of the 21 st 

278 century', while Cl y +aBr y is still above 1960 levels. Also, models with larger values of Cl y +aBr y 

279 (Figure 2c) do not necessarily have larger ozone loss, e.g., MR1 and CNRM-ACM have the 

280 largest peaks in ozone loss, but do not have the largest changes in Cl y +aBr y . 

281 One possible cause for the lack of a simple relationship between ozone and halogen 

282 changes is differences in temperature trends. As shown in Figure 2e, all models show cooling 

283 throughout the 21 st century, largely from increasing CO 2 concentrations, and this is expected to 

284 increase ozone (by slowing down the rate of reactions that destroy ozone). However, it should 

285 be noted that over the recent past about 50% of the cooling is associated with the ozone loss from 

286 Cl y +aBr y [ Shepherd and Jonsson, 2008]. This cooling explains the increase in ozone above 

287 1960 levels before Cl y +aBr y returns to 1960 values. Most models cool by about 8-9 K below 

288 1960 levels by the end of century, but a few models have cooling exceeding 10 K (CCSRNIES, 

289 UMSL1MCAT, LMDZrepro, and CMAM). This is consistent with an expected increase in 

290 ozone due to cooling, and these four models have among the largest ozone increases by 2100, as 

291 shown in Figure 2b. 



13 


292 To quantify the relative contributions of the above mechanisms to ozone changes we now 

293 apply an MLR analysis to model outputs. As described in Section 2, the MLR analysis is used to 

294 quantify the contributions of changes in Cl y +aBr y , temperature, NO y , and HO x to the ozone 

295 changes. As previously discussed in Section 2, there are a number of caveats with this type of 

296 analysis and care should be taken in interpreting the results. 

297 The contributions from Cl y +aBr y , temperature, NO y , and HO x , are shown in Figures 

298 2d,f,h,j, respectively. Differences among models in the contribution from Cl y +aBr y are 

299 generally similar to relative variations in Cl y +aBr y (compare Figures 2d and 2c), with some 

300 exceptions. Most notably, MRI and CNRM-ACM have the largest ozone loss at this level from 

301 1960-2000 despite not having unusually large amounts of Cl y +aBr y . This is related to larger 

302 sensitivities to Cl y +aBr y , and we quantify this sensitivity below. 

303 Temperature trends play an important role in the return of ozone to above historical 

304 levels: the four models with the largest cooling trend (LMDZrepro, CCSRNIES, CMAM, and 

305 UMSLIMCAT, see Figure 2e) have the largest increases in ozone relative to 1960 levels (see 

306 Figure 2b). The differences in cooling among models explain most of the differences in 

307 projected ozone by 2100. Although, again, we note that it is not a simple one-to-one 

308 relationship. 

309 The contributions from NO y and HO x changes are significantly smaller in most cases. 

310 Figure 2h,j show the contributions of NO y and HO x to ozone changes. For NO y most models 

311 show losses of 0.2-0.3 ppm by 2100 with a few models showing somewhat larger loss, up to a 

3 12 loss of 0.7 ppm in AMTRAC3. While 5 hPa is typically one of the largest loss regions from NO y 

313 changes in volume mixing ratio, HO x has a much smaller impact at this level and typically is 

314 more important at altitudes above 5 hPa [. Jackman et al., 1986]. All models show a relatively 
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315 small contribution (less than 0. 1 ppm) from HO x except ULAQ, which has the largest HO x trend 

316 (see Figure 2i). In general, for the Alb greenhouse gas scenario considered in these CCMVal-2 

317 simulations, NO y and HO x contributions are much smaller than the contributions from halogen 

318 recovery and stratospheric cooling over the 21 st century. It is however important to note that this 

319 may not be the case for other GF1G scenarios [e.g., Oman et al., 2010]. 

320 The above analysis indicates that changes in Cl y +aBr y and temperature dominate tropical 

321 ozone evolution at 5 hPa, and differences in the simulated changes in these quantities explain 

322 most (but not all) of the variations among the models. To see how representative the changes at 

323 5 hPa are over a range of pressures, Figure 3 shows vertical profiles of changes between 2000 

324 and 2100 for 25°S-25°N. 

325 Figure 3a shows each models reference ozone profile for the year 2000 while changes by 

326 2100 appear in Figure 3b. Also, plotted on Figure 3a are Microwave Limb Sounder (MLS) 

327 satellite measurements (solid black curve) averaged from 2005-2009. The peak ozone increase is 

328 typically at 3 hPa (Figure 3b) with most models centered at around 1.5 ppm. Figure 3c shows 

329 that in the upper stratosphere most models have a Cl y +aBr y decrease of 2 ppb over this time 

330 period. CCSRNIES and SOCOL have larger than average decreases in Cl y +aBr y and 

331 AMTRAC3 has a smaller than average Cl y +aBr y change with a different vertical structure, 

332 which is due to how AMTRAC3 treats the breakdown of halogens. Again, we see that MRI and 

333 CNRM-ACM have the largest ozone changes due to Cl y +aBr y (Figure 3d) with AMTRAC3 

334 having the smallest increase over the upper stratosphere. Most models show about equal 

335 contributions from decreases in Cl y +aBr y and decreases in temperature to the ozone increases 

336 over the 21 st century. Observations from the MLS [Santee et al., 2008a] for HC1 (not shown) 
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337 suggest that while a few models are in good agreement with MLS in the upper stratosphere, 

338 many appear to be on the low end with only CCSRNIES on the high end of observations. 

339 Figure 3 shows that the key features seen at 5 hPa hold throughout the upper stratosphere 

340 (10 - 1 hPa). Models that have high (low) ozone changes at 5hPa have similar characteristics 

341 throughout the upper stratosphere. Also, differences in Cl y +aBr y can explain most of the 

342 variations among the models with additional contributions from temperature. Figure 3h shows 

343 the contribution from NO y to changes in ozone: it is typically much smaller, and opposite in sign, 

344 than changes due to Cl y +aBr y or temperature. In general, models show NO y -related ozone loss 

345 of 0.2 ppm around 3-10 hPa, with only AMTRAC3 exhibiting a larger loss. As presented in the 

346 Introduction, stratospheric cooling increases NO y loss [. Rosenfield and Douglass, 1 998] causing a 

347 smaller increase (and even some decreases, see Figure 3g) in NO y over the 21 st century than 

348 would be expected from N 2 O increases in the absence of climate change. The interpretation of 

349 changes in lower stratospheric ozone are complicated by the increasingly important impact of 

350 transport which is not explicitly included in the MLR analysis and will be discussed in more 

351 detail in Section 3.4. 

352 As discussed above there are not always one-to-one relationships between changes in 

353 Cl y +aBr y and temperature and their contribution to ozone changes. One possible reason for this 

354 is the different sensitivities of ozone to Cl y +aBr y and temperature. Figure 4 shows the vertical 

355 profiles of the sensitivity of tropical ozone (25°S-25°N) to Cl y +aBr y , temperature, and NO y 

356 calculated over the period 1960-2100. Overall there is good agreement between the sensitivities 

357 of the various models. All models show a peak CI y +aBr y sensitivity at 3 hPa, mostly around - 

358 0.35 ppm/ppb. However, two models stand out with much larger sensitivities, CNRM-ACM (- 

359 0.6 ppm/ppb) and MRI (-0.8 ppm/ppb) at 3 hPa. 
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360 Further insight into the sensitivity of modeled ozone to Cl y +aBr y can be gained by 

361 examining the abundance of CIO and BrO archived by the various CCMs. Here, we provide a 

362 snapshot of the analysis presented in Chapter 6 of SPARC CCMVal [2010]. Figure 5 shows 

363 profiles of the C10/Cl y ratio archived by four of the CCMs at 35°N for Sept 1993 and at 22°N for 

364 Feb 1996 (colored lines). The black lines represent the calculated ratio of C10/Cl y found using a 

365 photochemical steady state (PSS) box model (e.g., Canty et al., 2006 and references therein) 

366 constrained by fields of O 3 , H 2 O, CH 4 , Cl y , NO y , aerosol surface area, and other constituents 

367 archived by each CCM group. Results are shown for four models: the two in question (MRI and 

368 CNRM-ACM) and two that archive CIO values that are well explained by the PSS simulation 

369 (CMAM and WACCM). The numerical value on each panel is the “metric” assigned to each 

370 comparison, with a value of unity indicating perfect agreement with the PSS simulation. 

371 The profiles of CIO archived by CMAM and WACCM are in excellent agreement with 

372 the PSS simulation, indicating that implementation of the chemical mechanism within these 

373 models agrees well with the implementation of the mechanism within the PSS model. However, 

374 the MRI CIO profile is much higher (green dashed line, Figure 5) than the profile found using the 

375 PSS model. As explained in SPARC CCMVal [2010], the MRI chemical mechanism neglected 

376 to include the channel for the loss of CIO by the reaction of CIO+OH to yield HCI+O 2 . This 

377 leads to much higher values of CIO than are reported by other models, or observed in the 

378 atmosphere [Santee et al., 2008b], Even when we neglect this reaction channel in our PSS 

379 model (black dotted curves. Figure 5), we still can not account for the high values of CIO 

380 archived by the MRI group. We conclude the high sensitivity of 0 3 to Cl y +aBr y exhibited by the 

381 MRI model is due to the high concentrations of CIO in this model, due in part to the neglect of 

3 82 the CIO+OH— HCI+O 2 channel. 
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383 The reason for the high sensitivity of O3 to Cl y +aBr y within the CNRM-ACM model is 

384 not known. For most altitudes, CIO within CNRM-ACM is simulated well by the PSS model. 

385 CNRM-ACM does simulate much higher values of CIO just above the tropopause than are 

386 accounted for by the PSS simulation; within CNRM-ACM, it appears that chlorine activation is 

387 occurring for warmer conditions than suggested by the PSS comparison. The high CIO just 

388 above the tropopause can not explain the large sensitivity of 0 3 within this model to Cl y +aBr y 

389 shown in Figures 2d, 3d, and 4a. We have examined profiles of BrO within the CNRM-ACM 

390 and the PSS model (not shown); we do not believe, based on this comparison, that the high 

391 sensitivity of 0 3 to halogens exhibited by CNRM-ACM is due to the calculation of BrO within 

392 this CCM. It is important to note that the PSS comparisons were done for the REF-B1 

393 simulation (1960 to 2000, observed aerosols) whereas the bulk of this analysis is based on the 

394 REF-B2 simulation (I960 to 2100, background aerosols). Finally, our comparisons cannot 

395 reveal possible errors in the coding of a rate constant that would affect the rate of ozone loss by 

396 CIO or BrO, but not affect the abundance of CIO or BrO (i.e., the reaction of CIO+O). The 

397 reason for high sensitivity of O3 to Cl y +cxBr y within the CNRM-ACM model remains 

398 unexplained. 

399 We performed an uncertainty analysis of the sensitivities as in Oman et al. [2010] and 

400 find similar results (not shown). That is, in the mid to upper stratosphere there is relatively more 

401 uncertainty in NO y sensitivity compared to much smaller uncertainty associated with Cl y +aBr y 

402 and temperature. 

403 The evolution of upper stratospheric ozone and of the different contributing ozone loss 

404 mechanisms at midlatitudes is very similar to that in the tropics (not shown). There are again, 

405 reasonably balanced, equal contributions of Cl y +aBr y and temperature for the midlatitudes as we 
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406 found in the tropics. Again MRI and CNRM-ACM stand out with much larger contributions 

407 from Cl y +aBr y to ozone changes. 

408 

409 3.4 Lower Stratosphere 

410 

411 As shown in Figure 1 the evolution of ozone in the tropical lower stratosphere (p > 20 

412 liPa) differs from that in the upper portion. This can also be seen in the profiles shown in Figure 

413 3b. 

414 Although an MLR analysis has been applied to the lower stratosphere, interpretation of 

415 these results is complicated by the large role of transport changes and the tight coupling of 

416 upwelling and temperature. In the tropical lower stratosphere the loss that the MLR attributes to 

417 a temperature change is actually largely a response to increased upwelling, which acts to 

418 decrease ozone. This is consistent with Avallone and Prather [1996], who showed that changes 

419 in upwelling are dominant in controlling tropical lower stratospheric ozone. Temperatures can 

420 be reduced by both increased in upwelling and decreased ozone. Previous studies have shown 

421 that an increase in the tropical upwelling is the principal cause of a decrease in tropical ozone 

422 [Shepherd, 2008; Li et al., 2009] with an additional possibly smaller contribution from “reverse 

423 self-healing”. Reverse self-healing occurs as upper stratospheric ozone rises above historical 

424 levels causing less ultra-violet radiation to penetrate into the lower stratosphere resulting in the 

425 decreased production of ozone. 

426 We examine here the relationship between tropical upwelling and ozone in the CCMVal- 

427 2 models. Figure 6 compares the change in tropical ozone (25°S-25°N) at 50 hPa between 1960 

428 and 2100 to the corresponding change in the vertical residual velocity (w*) at 70 hpa, for the 12 
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429 models that provided w* fields. All models show an increase in tropical upwelling and a 

430 decrease in ozone over this time period. However, there is a large spread amongst the models in 

43 1 the change in both quantities. There is a fairly compact relationship between these two quantities 

432 and a linear fit (black line) nearly intersects the origin indicating that increases in upwelling are 

433 the dominant contributor to ozone loss at this level. Most models indicate increases in upwelling 

434 of 0.04-0.10 mm/s over this time period and ozone decreases of 0. 15-0.35 ppm. SOCOL appears 

435 to stand out from this group with significantly larger increases in upwelling and ozone decreases, 

436 which act in concert to produce the very large temperature reductions seen in Figure 3e. The 

437 impact of this larger change in upwelling in SOCOL can also clearly be seen in both total 

438 column ozone change (see Figure lb), and partial lower stratospheric ozone change (see Figure 

439 lh). In general, differences in the change in upwelling explain differences in lower stratosphere 

440 tropical ozone in the CCMVal-2 models. Also, this linear relationship exists when using the 

441 same level for ozone and w* (both at 50 and 70 hPa). 

442 In the midlatitude lower stratosphere transport is also thought to play an important role 

443 for ozone changes, which are generally positive over the 21 st century, however, the exact role is 

444 not quantified in the present analysis (see further discussion below). 

445 

446 4. Polar ozone 

447 

448 We now consider changes in polar ozone. Austin et al. [2010b] evaluate the CCMVal-2 

449 model simulations of the Antarctic ozone hole in detail, and describe both the processes that are 

450 well simulated and the areas that need additional work. More information on model biases, 
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451 especially those involving temperature thresholds for PSC formation, are provided in Austin et 

452 al. [2010b], 

453 We will first focus on annual average changes and then look at a more detailed analysis 

454 of spring when the largest ozone depletion is observed. Figure 7 shows the evolution of annual 

455 average southern and northern polar ozone for total column ozone (upper panels, 7a,b), upper 

456 partial column (middle, 7c, d) and lower partial column (lower 7e,f). Again, we show ground- 

457 based total column ozone measurements (solid black curves, updated from Fioletov et al., 

458 [2002]) with respect to 1964 in Figures 7a, b. 

459 Evolution of the upper partial columns is very similar to that of midlatitudes (see Figures 

460 ld,f). Flowever, the evolutions of the lower and total columns differ. All models show larger 

461 peak ozone loss around 2000 in the southern polar region than in the northern polar region. This 

462 is expected and seen in observations, as there is greater polar stratospheric cloud (PSC) 

463 formation in the southern polar region and a more stable polar vortex that lasts into spring each 

464 year. In general there are larger increases in ozone amounts by 2100 compared to 1 960 in the 

465 northern polar region. It is difficult to say quantitatively with the existing simulations how much 

466 of the differences seen by 2100 in the total and partial lower stratospheric column ozone are due 

467 to polar chemistry differences, circulation differences, or a combination of the two. Additional 

468 simulations with fixed low chlorine concentrations or fixed GFIG concentrations could help to 

469 quantify the relative impacts. 

470 The largest ozone losses occur in the spring in the lower portion of the stratospheric 

471 column so we focus on this time and region in more detail. Figure 8a shows the Antarctic (60- 

472 90°S) lower stratospheric column ozone for October. Peak ozone loss averages about 100 DU 

473 around 2000 with a range from models between 70 and 170 DU. We applied a linear regression 
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474 analysis to the lower stratospheric column ozone amounts to calculate the contribution from 

475 Cl y +aBr y using a=60. as appropriate for the polar lower stratosphere. 

476 Figure 8b shows the contribution of Cl y +aBr y and the residual (Figure 8c), which is much 

477 smaller and for most models typically less than ±20 DU. Cl y +aBr y dominates the partial column 

478 ozone until the end of the 21 st century when the residual is in most cases of comparable 

479 magnitude. The decrease in mean age of air (i.e. stronger circulation, see Figure 8e) by itself and 

480 in combination with upper stratospheric ozone returning to above historical levels acted to cause 

481 increased ozone change (i.e, positive residual), however this could be tempered by the cooling of 

482 the polar lower stratospheric temperatures (Figure 8f shows 50 hPa temperatures) which occurs 

483 in nearly all models. 

484 Figure 9a shows the Arctic (60-90°N) lower stratospheric column ozone for March. We 

485 also note that this area is typically bigger then the size of the polar vortex so this average also 

486 includes mixing effects (also to a lesser extent in the Southern Hemisphere during October). 

487 Peak ozone loss averages at 30 DU around 2000 with a range between 10 and 60 DU. Again we 

488 linearly regressed the Cl y +aBr y against the partial lower stratospheric column ozone. Figure 9b 

489 shows the contribution of Cl y +aBr y and the residual (Figure 9c), which in this case is generally 

490 larger with a consistent 30-40 DU residual by 2100 for most models. This seems to indicate that 

491 in the NH polar regions circulation changes cause a larger increase in ozone than seen in the SH 

492 by 2100. 

493 

494 5. Ozone Recovery 


495 
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496 The above differences in ozone return dates between regions are summarized in Figure 

497 10 which shows the evolution of the multi-model mean (MMM) total and partial column ozone, 

498 with respect to 1960 for high latitudes (solid curves) and midlatitudes (dashed curves). Total 

499 column ozone for the MMM (Figure 10a) shows a return to 1960 values by the 2030s in the NH 

500 in both mid and high latitudes. In the SH the total column ozone MMM exhibits a return to 1960 

501 values by the 2060s in the midlatitudes and by 2080 at high latitudes. Nearly all the differences 

502 in total column ozone evolution are from the lower portion of the column (500-20 hPa, see 

503 Figure 10c) with very similar changes in both hemispheres occurring in the upper portion (Figure 

504 10b). 

505 Comparing with the evolution of halogens (red curves; using 60-90°S at 50 hPa as a 

506 proxy representative value), we see a significant difference in the dates ozone and halogens 

507 return to historical levels (i.e., 1960 levels). The MMM Cly+aBry does not return to I960 

508 values by 2100. This is later than the return of ozone (for extratropical regions). The cause of 

509 the earlier return of ozone in midlatitudes (shown in Figure 16 in Austin et al., 2010a) in both 

510 hemispheres is likely from an increased circulation [Austin and Wilson , 2006; Waugh et al., 

511 2009; Li et al., 2009] along with increased upper stratospheric ozone from stratospheric cooling. 

512 There is a noticeable difference in the dates when northern and southern hemisphere 

513 ozone returns to 1960 values (with no such difference in Cly+aBry, not shown). This 

514 interhemispheric difference is likely due to interhemi spheric differences in transport. Some of 

515 the delay in the SH return is likely from asymmetries in polar ozone loss since there is significant 

516 polar ozone loss occurring in the SH in the early 21 st century when NH ozone is returning to 

517 1980 levels (not shown with reference to 1980 levels), causing less ozone to be transported into 

5 1 8 SH midlatitudes when the vortex breaks down [ Atkinson et al., 1989]. The impact of polar ozone 
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519 losses being transported into the midlatitudes likely play the dominant role in the later return to 

520 historical values seen in the SH midlatitudes during first half of the 21 st century with any 

521 asymmetries in the strengthening of the stratospheric circulation likely having a larger impact 

522 toward the later half of the 21 st century'. 

523 

524 6. Conclusions 

525 Simulations of the past and future were preformed using 14 chemistry-climate models for 

526 the CCMVal-2 activity. Austin et al. [2010a] have examined the evolution of total column ozone 

527 and compared the recovery of column ozone and Cl y in these models. Here we have contrasted 

528 the evolution in the upper and lower stratosphere, and examined the cause of the ozone changes 

529 (and differences among the models) particularly for the 21 st century. 

530 The simulations presented here and in Austin et al. [2010a] show that there is some 

531 general agreement in the ozone evolution among the models, with all showing column ozone 

532 decreasing from 1960 to around 2000, and then increasing over the first half of the 21 st century. 

533 Models also show that the column ozone evolution varies with latitude, especially in the latter 

534 half of the 21 st century. There are clearly some quantitative differences in ozone evolution 

535 across the CCMVal-2 simulations, e.g., a large spread in simulated return of ozone to historical 

536 values [Austin et al, 2010a], 

537 Separation into partial columns above and below 20 hPa reveals that the latitudinal 

538 differences in the evolution of the column ozone are almost completely due to differences in the 

539 lower stratosphere (the region of the atmosphere below 20 hPa). In all models, there are only 

540 weak latitudinal variations in the evolution of upper stratospheric ozone, and at all latitudes 

541 upper stratospheric ozone increases throughout the 21 st century and returns to 1960 levels before 
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542 the end of the century. There is, however, a large spread in dates of return to historical values, 

543 with, for example, dates of return to 1960 levels at 5 hPa varying from 2020s to 2060s. A 

544 multiple linear regression analysis indicates that the upper stratospheric ozone increase comes 

545 from almost equal contributions of decrease in halogens (Cl y +cxBr y ) and cooling from increased 

546 greenhouse gas concentrations (for the Alb greenhouse gas scenario considered in these 

547 simulations), with only small contributions from NO y and HO x . 

548 In the tropical lower stratosphere there is a steady decrease in ozone through the 21 st 

549 century in all models, whereas ozone in middle and high latitudes increases during the 21 st 

550 century. As a consequence tropical column ozone does not return to 1960 levels, whereas 

551 extratropical ozone returns or is close to 1960 levels before the end of the century. The decrease 

552 of tropical lower stratospheric ozone is due to an increase in tropical upwelling. Models with 

553 larger ozone decreases have larger increases in upwelling. Changes in transport also play an 

554 important role in the evolution of mid-latitude lower stratospheric ozone, and strongly contribute 

555 to the earlier return of ozone than Cl y +aBr y to historical values. 

556 There are quantitative differences between the hemispheres in the simulated ozone 

557 evolution. In particular, ozone returns to historical values earlier in the northern hemisphere. This 

558 difference in hemispheric return dates varies between 10 and 50 years depending on the latitude 

559 of interest and the reference year (e.g. 1960 or 1980) chosen, and is almost completely related to 

560 differences in the lower stratosphere. Additional simulations would need to be completed in 

561 order to quantitatively attribute causes for hemispheric differences. However, based on existing 

562 simulations it appears that the larger ozone loss from polar chemistry in the SH and the transport 

563 of this ozone depleted air into the midlatitudes as the vortex breaks down is one cause, especially 
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564 during the first half of the 21 st century. Changes in the stratospheric circulation could also play a 

565 role in these differences. 

566 At high latitudes, especially in the SH (Figure 8) Cl y +aBr y dominates long-term trends in 

567 the lower portion of the column ozone (and total column ozone) with a small (generally positive) 

568 residual. This residual is likely caused by circulation changes and upper stratospheric return to 

569 above historical levels. 

570 These results reinforce the conclusions in Eyring et al. [2007] using simulations from the 

571 CCMVal-1 activity: decreasing levels of halogens, continued stratospheric cooling, and changes 

572 in circulation are the major factors driving 21 st century ozone trends. It is important to keep in 

573 mind that all of these simulations are based on a single GFIG and halogen scenario, with only one 

574 model including an interactive ocean. In addition, the use of mixing ratio based halogen 

575 boundary conditions rather than emission-based constrains the models response [Douglass et al, 

576 2008]. Future work should be done to explore these issues and their impact on ozone evolution 

577 in more detail. 

578 
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851 Table 1. Model description and references. 


Model 

Atmospheric 

GCM 

Domain/Resolution 
or Truncation 

SST/Sea Ice for 
Ref-B2 

Model Reference 

AMTRAC3 

AM3 

Variable, ~200 km, 
48 L, 0.017 hPa 

CM2.1 

Austin and Wilson 
(2010c) 

CAM3.5 

CAM 

1.9° x 2.5°, 26 L, 3.5 
hPa 

CCSM3 

Lamarque et al. 
(2008) 

CCSRNIES 

Hlil 

T42, 34 L, 0.012 hPa 

MIROC / IPCC- 
AR4 

Akiyoshi et al 
(2009) 

CMAM 

AGCM3 

T31, 71 L, 0.00081 
hPa 

Interactive 

■BH 

CNRM-ACM 

ARPEGE- 
Climate 
version 4.6 

T42, 60 L, 0.07 hPa 

CNRM-CM3 

AR4 

Deque (2007); 
Teyssedre et al. 
(2007) 

GEOSCCM 

GEOS 5 

2° x 2.5°, 72 L, 
0.015 hPa 

HadlSSTl for 
Ref-Bl and 
CCSM3 for 
Ref-B2 

Pawson et al. 
(2008) 

LMDZrepro 

LMDZ 

2.5° x 3.75°, 50 L, 
0.07 hPa 

OPA (ocean), 
LIM (ice) 

Jourdain et al. 
(2008) 

MRI 

MJ98 

T42, 68 L, 0.01 hPa 

MRI- 

CGCM2.3.2 

Shibata and Deushi 
(2008a, b) 

SOCOL 

MAECHAM4 

T30, 39 L, 0.01 hPa 

ECHAM5- 

MPIOM 

Schraner et al. 
(2008) 

ULAQ 

ULAQ-GCM 

R6/1 1.5° x 22.5°, 26 
L, 0.04 hPa 

CCSM3 

Pitari et al. (2002) 

UMSLIMCAT 

HadAM3 L64 

2.5° x 3.75°, 64 L, 
0.01 hPa 

HadGEMl 

Tian and 
Chipperfield 
(2005); Tian et al. 
(2006) 

UMUKCA- 

METO 

HadGEM-A 

2.5° x 3.75°, 60 L, 
84 km 

HadGEMl 

Davies et al. 
(2005); Martin et 
al. (2006); 
Morgenstem et al. 
(2009) 

UMUKCA- 

UCAM 

HadGEM-A 

2.5° x 3.75°, 60 L, 
84 km 

HadGEMl 

Davies et al. 
(2005); Martin et 
al. (2006); 
Morgenstem et al. 
(2009) 

WACOM 

CAM 

1.9° x 2.5°, 66 L, 
0.00000596 hPa 

CCSM3 

Garcia et al. (2007) 


852 


Figures 



























34 



854 Figure 1. Annual average total and partial column ozone amounts over 3 regions including the 

855 tropics (25 0 S-25°N) and midlatitudes of each hemisphere (35-60°S and °N). The partial column 

856 ozone is separated into an upper portion (d,e,f) from 20 to 0. 1 hPa and a lower portion (g,h,i) 

857 from 500 to 20 hPa. All data are from 1960 to 2100 (except UMUKCA-METO to 2083) and 

858 have been smoothed with a 1:2:1 filter iteratively 30 times. Note the scale change in y-axis in 

859 panels d,e,f which are half the magnitude of the other panels. Ground-based total column ozone 

860 observations (black curves) are shown over 1964-2007 relative to 1964. 
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Figure 2. Evolution of ozone (a) in the tropics (25°S-25°N) at 5 hPa and change in ozone (b) 
with respect to 1960 levels for the CCMVal models. Models change in Cl y +5Br y (c) and 
temperature (e) with respect to 1960 levels. Contribution of changes in Cl y +5Br y (d) and 
temperature (f) to changes in ozone. All models are shown from 1960 to 2100 (except 
UMUKCA-METO to 2083) and have been smoothed with a 1 :2:1 filter iteratively 30 times. 
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Figure 2 (cont.). Model evolution of NO y (g) and HO x (i) with respect to 1960 levels over (25°S- 
25°N) at 5 hPa. Contribution of changes in NO y (h) and HO x (j) to changes in ozone. All 
models are shown from 1960 to 2100 (except UMUKCA-METO to 2083) and ha\e been 
smoothed with a 1:2:1 filter iteratively 30 times. 
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Figure 3. Profiles of ozone (a) in the tropics (25°S-25°N) for 2000 and differences in ozone (b) 
from 2000-2100 for the CCMVal models. Model differences in Cl y +5Br y (c), temperature (e), 
and NO y (g) from 2000 to 2100. Also shown are the contributions of changes in Cl y +5Br y (d), 
temperature (f), and NO y (h) to changes in ozone. Except UMUKCA-METO change shown 
from 2000 to 2083. Also shown in (a) are observations from MLS (solid black curve) for 
reference. 
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Figure 4. Tropical (25°S-25°N) profiles of the sensitivity of ozone to changes in (a) Cl y +5Br y 
(ppm/ppb), (b) temperature (ppm/K), and (c) NO y (ppm/ppb). 
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890 

891 Figure 5. Comparison of zonal, monthly mean profiles of radicals from CCM models (coloured 

892 lines) versus 24-hour average radical profiles found using a PSS box model constrained by 

893 profiles of T, O 3 , H 2 O, CH 4 , CO, NOy, Cly, Bry, and sulfate SAD from the various CCMs for 

894 the two indicated times and locations. The PSS model was run at CCM model levels from the 

895 tropopause (dashed colored lines) to 1 hPa. All simulations used JPL 2006 kinetics. The colored 

896 error bars represent the standard deviation about the zonal monthly mean for various days used 

897 to compute the mean. The black error bars represent the sensitivity of PSS output to variability in 

898 the CCM profiles of radical precursors. The metrics (“values of g”) were found as described in 

899 Chapter 6 of SPARC CCMVal2 [2010], For the MRl model. PSS simulations are shown with 

900 and without loss of CIO by the reaction CIO+OH— MICI+O 2 . 

901 
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70 hPa Aw* vs. 50 hPa AO. (25S-25N, Ann.) 1960-2100 
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907 Figure 6. Scatter plot showing the change (from 1960-2100) in 70 hPa w* and 50 hPa ozone. 

908 The values are annual averages over 25°S-25°N for 12 of the CCMVal models with the black 

909 line showing the linear fit. UMUKCA-METO change shown from 1960-2083. 

910 
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912 Figure 7. Annual average total and partial column ozone amounts over the high latitudes of each 

913 hemisphere (60-90°S and °N). The partial column ozone amounts are separated into an upper 

914 portion (c,d) from 20 to 0.1 hPa and a lower portion (e,f) from 500 to 20 hPa. All shown from 

915 1960 to 2100 (except UMUKCA-METO to 2083) and have been smoothed with a 1:2:1 fdter 

916 iteratively 30 times. Note the scale change in y-axis in c,d. Ground-based total column ozone 

917 observations (black curves) are shown from 1964-2007 with respect to 1964. 

918 
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Figure 8. Evolution of ozone (a) for October partial column ozone (500-20 hPa) over 60°-90°S 
with respect to 1960 levels for the CCMVal models. Model contribution of Cl y +60Br y at 50 hPa 
(b) and residual (c) with respect to 1960 levels. Changes in Cl y +60Br y at 50 hPa (d), change in 
mean age of air (e) for models that included age of air tracer, and (e) change in temperature. All 
shown from 1960 to 2100 (except UMUKCA-METO to 2083) and have been smoothed with a 
1:2:1 filter iteratively 30 times. 
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Figure 9. Evolution of ozone (a) for March partial column ozone (500-20 hPa) over 60°-90°N 
with respect to 1960 levels for the CCMVal models. Model contribution of Cl y +60Br y at 50 hPa 
(b) and residual (c) with respect to 1960 levels. Changes in Cl y +60Br y at 50 hPa (d), change in 
mean age of air (e) for models that included age of air tracer, and change in temperature (f). All 
shown from 1960 to 2100 (except UMUKCA-METO to 2083) and have been smoothed with a 
1:2:1 filter iteratively 30 times. 
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Figure 10. Multi-model mean (MMM) of total (a) and partial column (b,c) ozone amounts for 
mid (35-60°S and 35-60°N, dashed curves) and high latitudes (60-90°S and 60-90°N, solid 
curves) and Cly+aBry for high latitudes (60-90°S) for reference. 






